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Abstract 

We review a perturbative approach to deal with Lagrangians with higher or 
infinite order time derivatives. It enables us to construct a consistent Poisson 
structure and Hamiltonian with only first time derivatives order by order in 
coupling. We show that, to the lowest order, the Hamiltonian is bounded from 
below whenever the potential is. We consider spacetime noncommutative field 
theory as an example. 



1 Introduction 



In the canonical formulation, we usually consider Lagrangians with only first time 
derivatives. However, higher derivative theories, including nonlocal theories, also have 
many physical applications. For example, when one integrates out high energy degrees 
of freedom in a local field theory, the low energy effective action is generically nonlocal 
[0. Higher derivative theories were considered even earlier in order to find a finite 
quantum field theory 0, before the advent of renormalization. Moreover, theories with 
infinite derivatives are inevitable from the viewpoint of string theory 0, ^. There are 
many other examples, such as higher derivative gravity , meson-nucleon interactions 
0, and spacetime noncommutative field theory |]^, |[, and so on. 

Unfortunately, an exact treatment of generic nonlocal theories is so far unknown. 
This was viewed as a serious problem for string theory However, in most cases, 
higher derivative terms appear as higher order corrections in the effective Lagrangian, 
so a perturbative approximation scheme would already be very useful. 

In this paper we consider theories for which the free part of the Lagrangian involves 
only first time derivatives, but there can be higher derivatives to any (finite or infinite) 
order in the interaction terms. According to the canonical formulation, the phase space 
is 2k dimensional for Lagrangians with fc-th time derivatives 0. For a low energy 
effective theory, however, we usually assume that there are as many degrees of freedom 
as the free part of the Lagrangian tells us. This is also an assumption for the S-matrix 
to be unitary in a quantum field theory |]lO|. Presumably the extra dimensions of 
the phase space are associated with certain high energy degrees of freedom which are 
integrated out from the description. 

What we need is a formalism to deal with Lagrangians with finite or infinite deriva- 
tives without introducing new (high energy) degrees of freedom. It should also allow 
us to consistently construct Poisson brackets and Hamiltonian perturbatively order by 
order, and in the end give us an effective Lagrangian with only first time derivatives. 



In |Tl]], it was proposed to remove higher derivative terms in the Lagrangian by field 
redefinition. But it is unclear how to achieve this to arbitrary orders for a generic theory, 
or whether this is always possible. A general procedure which provides a Hamiltonian 
formulation order by order in the number of time derivatives in the interaction terms 
was constructed in |jl2|. The perturbative approach [Q that we will discuss in this paper 
is a modification of that. It is based on an expansion in the coupling constant, with no 
restriction to the number of derivatives in each order. This technique has been used 



extensively in [|T^ for higher derivative gravitational theory. Although this approach 
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was studied in detail in which we recommend the interested readers to refer to for 
related issues, such as the quantization or causality problems, we present this approach 
here in a slightly different fashion, with some derivations and proof explicitly given. In 
fact, we have discovered the same techniques independently, and were only acquainted 
with previous works on this subject by a referee after this work was finished. 

What is the new progress we made on this subect? (1) We show that, to the lowest 
order, we can describe the result most simply by a change of variable. Remarkably, 
it turns out that, in terms of this new variable, the Hamiltonian is bounded from 
below whenever the potential is. (2) We give the formal proof that the perturbative 
approach can be carried out to all orders consistently. This proof also shows that if an 
exact solution to the equation of motion is known, the same approach can immediately 
give the final result. (3) In order to understand the difference between the original 
system and its perturbative formulation, we demonstrate by an example that when 
the calculation is carried out to all orders, the description of low energy modes is 
exact, including the effect of their interactions with high energy modes, but high energy 
degrees of freedom are still absent. (4) We apply this approach to field theories living 
on noncommutative spacetime, and comment on the spacetime uncertainty relation. 

This paper is organized as follows. In Sec. 0, we briefly review the canonical 
formulation of theories with finite time derivatives, and explain the problems in dealing 
with infinite derivatives. We comment there and show in the appendix that one-loop 
unitarity is not affected by introducing higher spatial derivatives, but is at stake when 
higher time derivatives are present. In Sec. |^ the perturbative approach for higher 
derivative theories is constructed. One can have a consistent Hamiltonian formulation 
of the theory to an arbitrary order in coupling. We also give a formal proof that this 
procedure can be consistently carried out to all orders. Finally, in Sec. ^ we give three 
examples. In the first example 0] we demonstrate explicitly how the perturbative 
approach works to the second order. The second example is used to show that this 
approach suppresses high energy degrees of freedom, but provides an exact description 
for the low energy modes. We also apply this approach to the example of spacetime 
noncommutative field theory. 



2 Canonical Formulation 
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2.1 Lagrangians with Finite Time Derivatives 



Here we review the canonical formulation for classical theories with finite time deriva- 
tives. Consider a Lagrangian with fc-th time derivatives 

Lo = Lo(g,g,g, (1) 

where g'-*-' is the i-th time derivative of q. Let us apply the formalism of The 
variation of the action 



S 



dt L 



with respect to q is found to be the time integral of 5q times the Euler-Lagrange 
equation 



±4. 



dLn 



n=0 



plus a boundary term 



dt I (9g(") 



0, 



(2) 



(3) 



where Pi is the conjugate momentum of g^*^ 



Pi 



k-i-l 

E 

j=0 



d 



dt) 9g(*-^+i)' 

The symplectic structure can be directly read off from this as 



(4) 



k-l 

n = Y^dPidq^^. 

i=0 



(5) 



The Hamiltonian is 



k-l 



H = J2P,q^^^-L{q,q,---,q^'^-'\q 



(6) 



2=0 



where g^^-* is viewed as a function of g'^^^'s for i = 0, ■ ■ ■ , k — 1 and Pk-i- Q Since H is 

linear in Pi, ■ ■ ■ , Pk-2, the energy spectrum is unbounded from below. 

Another approach to deal with higher derivative Lagrangians is to introduce new 

variables to replace higher derivatives of g, and use Lagrange multipliers to ensure 

that we have not introduced any new degrees of freedom. One natural choice is to 

define new variables g^ by g^ = g*^*^ (Exercise 1.26 in [jl5l), and another is to set g^ = 

^ We assume that the Lagrangian is nondegenerate. This can always be achieved by adding total 
derivatives. 
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g*^^*). In any case, after introducing new variables, the Lagrangian only has first time 
derivatives. Canonical momenta of all variables, including the Lagrange multipliers, can 
be defined as usual. Some of these will impose primary constraints. Dirac's procedure 
for constrained quantization can then be used to obtain the same result. See |jl5[ for a 
detailed description of the general procedure. 

2.2 Lagrangians with Infinite Time Derivatives 

For a Lagrangian with infinite time derivatives Lo{q, q,- ■ ■), one may hope that a trick 
similar to the one for Lagrangians with finite derivatives will work. A natural guess is 
to define the new Lagrangian as [ITU 



L = LoiQ) + / dX /i(A)(g(A) - g'(A)), (7) 



where Lq{Q) is obtained from Lo[g] by the replacement 

q{t)^Q{t,0), gW(t)^_Q(t,A)U=o. (8) 

By Q' we mean -^Q and is a Lagrange multiplier. This implies that Q{t, A) = qit + X). 

However, the Euler-Lagrange equations for this new Lagrangian L is in fact different 
from the Eular-Lagrange equation for the original Lagrangian Lq. Further constraints, 
or boundary conditions, have to be imposed |l6|-|jl9 . 



Instead of trying to resolve this problem, we give some general reasons why a naive 
extension of the idea for finite time derivatives has difficulties. 

For a nondengerate Lagrangian with k-th time derivatives, the phase space is 2k 
dimensional for a single variable q. Let us choose g*^*-* (z = 0, 1, ■ ■ ■ , 2A; — 1) to be the 
coordinates on the phase space. The symplectic 2-form (^) implies that the Poisson 
bracket 

(g«,g(^))=0 (9) 

for any i,j < k — 1. The bracket is nontrivial when i < k — 1 and j > k and vise versa. 
When we take k ^ oo, naively, the phase space has coordinates {q^^^i=o,i, --,oo, and 
holds for any i and j. The nontrivial part of Poisson brackets is now completely 

lost! 

On the other hand, if we use g^*^ and Pi (i = 0, 1, ■ ■ ■ , — 1) as the phase space 
variables, in the limit k ^ oo, the phase space has coordinates -Pi}i=o,i,---,oo, which 
seems to have twice the dimension of the previous guess. 

Another problem with infinite derivatives is the Cauchy problem |2^. For a La- 



grangian with k-th derivative its equation of motion is a differential equation of order 
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2k. To determine the trajectory uniquely, the initial values of g*-*-* for i = 0, 1, ■ ■ ■ , 2/;; — 1, 
have to be given. Naively, when k oo, all g^^^'s have to be given as initial conditions, 
but they will completely determine the trajectory of q via Taylor expansion even with- 
out using the equation of motion. This seems to suggest that for infinite k the Euler 
Lagrange equation should be imposed as a constraint on the space of initial conditions, 
which is equivalent to the phase space, and the evolution for a point in the phase space 
is trivial. This viewpoint is consistent with the vanishing of Poisson brackets mentioned 



above, and is similar to the proposal of earlier works |1I6|-|[T9[ . 

Although it is desirable to have a general formalism to deal with theories with 
infinite derivatives. It is known that different cases can have very different properties. 



For some interesting examples see 0,|17|- Probably different theories with infinite 



derivatives may have to be treated very differently. 

Normally, one needs to understand the classical Hamiltonian formulation of a system 
before quantization can be done. In the operator formalism we replace the Poisson 
brackets by commutators. In the path integral formalism, we should start with 

J Dq Dp el '^'^P'^-"^ (10) 

for a system with unconstrained conjugate variables p, q. The formula 

l^ge/''*^ (11) 

may or may not be correct. 

Given a Lagrangian with infinite derivatives, if one simply use (|Tl|) to define a 
quantum system^, the violation of unitarity is observed in several examples |21]. It is 



tempting to interpret Q(t, A) in (|^ as a chiral field on an open string, and an intriguing 
possibility is that for certain inifinite derivative theories one can find a procedure to 
add new degrees of freedom for consistency which eventually leads to string theory. 

In the next section we will review a perturbative approach for low energy effec- 
tive theories with infinite time derivatives. This approach allows us to change the 
Lagrangian order by order to a new Lagrangian with only first time derivatives. It 
is then possible to quantize the theory without breaking unitarity. In the appendix 
we also show that, with the free part of the Lagrangian unchanged, arbitrary spatial 
derivatives in the interaction terms will not spoil unitarity at one-loop level, which is 
sufficient for low energy effective theories. In particular, an immediate consequence is 
that noncommutative field theories with only spatial noncommutativity is unitary at 
one loop level. 

^ The Hamiltonian formulation of IT§1 also leads to (IT 
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3 Pert ur bat ive Approach 



In this section we review an approach to deal with higher derivative terms perturba- 
tively. This approach works equally well for Lagrangians with finite or infinite order 
time derivatives. 

In many applications, higher derivative terms appear in the Lagrangian as higher 
order corrections to a low energy effective theory. However, according to the canonical 
Hamiltonian formulation, when higher order terms arc introduced, the dimension of 
the phase space increases, no matter how small the correction term is. Intuitively, 
this discrete jump in dimension does not seem to be compatible with the physical 
interpretation of the higher derivative interaction terms as a perturbative correction to 
the theory. 

To resolve this apparent conceptual conflict, we note that the new dimensions of 
the phase space correspond to those degrees of freedom which are not accessible at 
low energies. If all we want is to have some knowledge about how these high energy 
states affect the behavior of the low energy modes, we can choose to restrict ourselves 
to the low energy degrees of freedom on the phase space. Therefore, the main idea 
behind this approach is to project the symplectic structure on the whole phase space 
onto the subspace of low energy states, and then induce the Hamiltonian formulation 
for the reduced phase space. This approach can also be generalized to cases in which it 
is preferred to work on a bigger phase space with finitely many higher time derivatives 
kept in the end. 

It is sufficient to illustrate this approach by considering a + 1 dimensional theory 
with a general Lagrangian 



To be definite, let us consider the cases with infinite order time derivatives. Under 
variation, the action is 






(13) 



where the equation of motion is 




(14) 



and Pn is the canonical momenta for q^'^^ 




(15) 



6 



The symplectic two-form is 

oo 

Q = Y.dPkdq^''\ (16) 
3.1 First Order Approximation 

In an exact treatment, all g'^^^'s are independent. In the low energy approximation, we 
only keep q and q. Using the equation of motion, we can replace higher derivatives of 
q by lower derivatives. To the lowest order approximation, 



n = even 



^ (_^2)(n-l)/2- (^ = odd). ^ ^ 



The boundary term of 6S (|T^) reduces to 

[UoSq + U^Sq]lf = [{q - g^o)5q - g^M^ , (18) 

where 

oo oo / , \ n 2k 2 ot^ 

- S S (20) 



A;=0n=2fc+2 



Note that and should be viewed as functions of q and q via the replacement (|T^. 
We find the symplectic two-form ( p!6|) as 

(21) 



= dUodq + dUidq = ( — 1 + g-^ — g-^] dqdq. 

\ oq oq J 



It follows that the Poisson bracket should be 

{q,q)^l+g^-g^. (22) 
To the lowest order in g, the Poisson bracket can be written as 

{x,p) = 1, (23) 

where 

x = q + g^i, p = q- g^o- (24) 

The Hamiltonian is 

H = Uoq + n,q-L, (25) 
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where we replace all higher derivatives of q by (]T7|). We find 



1 Tfl ~ 

H=-p^ + —x' + gVix,p), (26) 

where V{x,p) = V{x,p, —m'^x, ■ ■ ■). Note that, if the potential V is bounded from 
below, the first order Hamiltonian is also bounded from below. 
One can check that the Hamilton equations 

q = {q,H), q={q,H) (27) 

reproduce the equation of motion (IHp to the first order in g. The corresponding 
Lagrangian is 

L = px — H = -X —X —gV{x^xj. (28) 

Its Euler-Lagrange equation agrees with the original system to the first order in g. Since 
the final expression only contains first derivatives, its quantization is straightforward. 

3.2 Higher Order Approximation 



For higher order corrections, we first iterate the equation of motion (|Tj) up to a certain 
order 0{g"'). For example, at the first order. 



where we should also use the replacement (|1^) to change the last term into a function 
of q and q only. Higher derivatives of q can also be replaced by functions of q and q 
up to the same order in g by differentiating with respect to time and repeatedly using 
m as 



H ^ I {-m^T'S -gY.tl ET=oi-mr^'^^ (1)'^'"' ^ (n = even), 

' ^[i-rn'Y-''/'<i-9Et,'^^'ET=oi-rn'Y^^^^^^ (n = odd). 

(30) 

In general, we can have all q^^^^s expressed as a function of q and q only, up to a 
certain order 0{g^). This helps us to derived the symplectic form from (p6[) 



— + ^ dqdq. (31) 
dq dq J 

The final Hamiltonian is defined by (|25|) with all higher derivatives of q replaced 
by functions of q, q. Its Hamilton equations (^71) gives the equation of motion up to 
0{g^^'). 
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An effective Lagrangian can be defined as 

L=px-H (32) 

for a conjugate pair {x,p) = 1. For sucli Lagrangians we do not expect violation of 
unitarity since it contains only first time derivatives. Although we expect such theories 
to be nonrenormalizable in the usual (power-counting) sense, it remains to be seen 
whether it is renormalizable in the mordern sense p^ . 



3.3 To All Orders: A Formal Proof 

Now we give a formal proof for the self-consistency of the perturbative formulation. 
From the equation of motion (p^, in principle one can find 



Q = f{<lA) 



(33) 



for a certain function / to all orders in g by iteration (or inspiration). Higher derivatives 
of q can be derived from this 

9^"^ = /n(g,g), (34) 

where the functions /„ can be obtained recursively 



n+l 



-f 



dfn . dfn . 

oq oq 



where we used the notation 



(35) 



(36) 



Here are a few identities that will come in handy in the following. From (|TBp, we 



find 



(37) 



For an arbitrary function A on the phase space, we also find the following identities 



1m 

dq 



fe=0 



d\A] 
dq 
OA 



dfk 
dq 



(38) 
(39) 



where one can also replace ^ by ^ in the last formula. 
The variables Hq, Hi can be read off from 



Uo6q + U^6q = Y.Pk6q^^\ 

k=0 



(40) 
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and we find 



Hi 



E 

A:=0 
oo 

E 

A;=0 L 



Pi 



Pi 



dq 
dq 



The Hamiltonian is 



H =[Uoq + U,q-L]. 



The Hamilton equations ([27|) based on the symplectic structure ( pTD are 



no + ni|f-f^(g-/) 



0, 



ni + no + n,|-f^(g-/)-^ 



0. 



(41) 
(42) 

(43) 

(44) 
(45) 



With the help of (|37| ) -(|39 D , one can show from (^11) , (|4^) that (^) is automatically 
satisfied, and that ( jl^ ) is equivalent to the equation of motion (^5]). 

Certainly, (|33| ) is not the only solution of the equation of motion (p^). It may be 
extended to a more general solution with many free parameters 



Q = f{Q,Q,Cl,C2,■■■)■ 



m 



These extra parameters are associated with high energy degrees of freedom. Ignoring 
them, or setting them to zeros corresponds to taking the low energy limit in which 
these high energy modes are "minimally excited". That is, they are only excited to 
be consistent with the equation of motion, but we do not consider their independent 
degrees of freedom. We will consider a system of coupled springs in the next section as 
an example to show this property more concretely. 



4 Examples 

4.1 An Example of Second Order Approximation 

In this section we show the perturbative procedure explicitly for an example to the 
second order. Consider the Lagrangian 

^ = + + gqf- (47) 

Its Euler-Lagrange equation is 

q = aq + gf + 2g{qqy. (48) 
10 



Following the procedure outlined in the previous section, we find 

no = g-2(7(gg)-, Tl^ = 2gqq (49) 
from (0), (0). Iterating p8|) to the first order in (7, we find 

g = ag + ^(5aV + W) (50) 

Taking derivatives, we get 

g(3) = aq + g{l^a''qq)+0{g''), 

= a2g + ^(23aV + 22aV) + 0(/), (51) 

and so on. 

Up to third order terms in we find 

Ho = g-4(7(agg) -2/(23aVg + 4ag3), (52) 
Hi = 2c/(ag2) +2/(5aV + 4agg^). (53) 

The Poisson bracket is 

(g, g) = 1 + ^gaq + 4^2(35^2^2 ^ g^^2^) ^ (>^(^3^_ (54^ 

The Hamiltonian is 

if = - + (/(a^g^ - 4agg2) + 2/(5aV - lOa^g^g^ - Aaq^) + 0(^^). (55) 



One can check that the Hamilton equations (|2^) give the equation of motion (|48|) correct 
to the third order in (7. 

The effective Lagrangian is 

L = pq-H, (56) 

where p is the phase space variable conjugate to q, i.e., (g,p) = 1. From (^), we can 
solve for p 

p = g - 8(7agg - 4g\l9a''q''q + ^aq^) + 0{g^). (57) 
Thus we find the Lagrangian for q, q 

L = -q^ + -q^ - g{a\^ + Aaqf) + 2g\-5a^q^ - 19a^qY - -aq^) + 0{g^). (58) 
2 2 3 

The Euler-Lagrange equation derived from this Lagrangian is 

q = aq + g{ha^q^ + Aaq^) + 4^2(19^3^3 ^ ^^a^qq^) + o{g^). (59) 

This is exactly what one gets from (^) by iteration, in the same way we got (0), but 
to a higher order correction in g. 

11 



4.2 Coupled Springs 

Consider two springs A and B of spring constants k and K, respectively, as in Fig. 1. 
One end of A is fixed to a wall. The other end of A is connected to B, with a mass M 
attached to the joint. Another mass m is attached to the end of B. 



Figure 1 

For K ^ k and M m, the degrees of freedom of M is much harder to excite than 
m. Let the coordinates of m and M he x and y. The equations of motion are 

mx + k{x — y) = 0, (60) 

My + Ky + k{y - x) = 0. (61) 

They imply that 

ax^^) + bx + cx = 0, (62) 

where a = Mm, b = (M + m)k + rriK, c = Kk. The effective Lagrangian for this 
equation of motion is 

L oc -(— ax^ + bx'^ — cx^). (63) 

Assume that K ^ k, then the first term is negligible for small x. Imagine that some- 
one studies the behavior of the spring system at low energies by conducting experiments 
in which only the variable x is manipulated, and finds an approximate Lagrangian as 
( |63D without the first term. The prediction of the natural frequence is 

. Kk k ( k M + m\ 

° {M + m)k + mK m\ K m ) ^ ' 

to the lowest order in 0{k/K). 

Now suppose that via more accurate experiments, the first term in the Lagrangian 
( |63[ ) is put in. When the perturbative approach is carried out to the first order, one 
finds 
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This is a better approximation than ( [6^ ) to the real natural frequency. In this case we 
can even ask what will we get by carrying out the perturbative approach to all orders. 
It is straightforward to see that the answer is 

= -4ac). (66) 

2a 

This is precisely one of the two normal mode frequencies of the exact equation of 
motion. Although it does not tell us anything about the other natrual frequency at 
higher energy u;^ ~ ^, it describes the exact effect of the stiff spring on the low energy 
mode. 

In this example, we see that the physical meaning of reducing the phase space to 
the space of g, q in the perturbative formulation is to suppress the excitation of high 
energy modes in a way consistent with the exact evolution, and only to allow low energy 
modes to be excited. 



4.3 SpaceTime NonCommutative Field Theory 

Spacetime noncommutative field theories are interesting and important for many rea- 
sons . Although the naive quantization is non- unitary and non-causal , perhaps 
a proper treatment will lead to the noncommutative open string (NCOS) theory [p3| , |2^. 

The exact treatment of this system is outside the scope of this paper. We will only 
study this theory perturbatively. The zeroth order results are in fact already studied 
in I^S], where it was shown that the kinetic term with or without the *-product results 
in the same Hamiltonian formulation. 

Consider the Lagrangian 

L = ^(9^09^0 -mV) + A0* 0*0. (67) 
The *-product is defined by 

f*g{x)=e¥'''^^^'^f{x)g{x')U,=,. (68) 

It makes no difference whether we put *-products in the quadratic terms. 
For simplicity, consider the 1+1 dimensional case 

[t,x], = e. (69) 

As we saw in Sec. p.l| , by a field redefinition, the effective Lagrangian is given by ( P5| ) 
at the lowest order. For the case at hand it is 

L ^ ^(a^$a^$ - m^^^) + A$*$*$, (70) 
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where the *-product is defined as * in ( |68D with the replacement 

f (-a2 + m2)"/2 (n = even), 

' ^{{-dl + m'Y^-^y^dt (n = odd). ^ ^ 

In the action defined by (^), one only needs 5"-terms for even n. 

This result clarifies a conceptual puzzle. On one hand, the uncertainty relation 

At Ax > 1^1, (72) 

which follows from (|69|) , implies that when we treat t as a classical number, Ax oo. 
Thus we should not associate independent physical degrees of freedom at different 



values of x. It suggests a reduction of the number of degrees of freedom p6|, |2^. On 
the other hand, higher derivatives in the action imply a larger phase space in canonical 
formulation. These two intuitions seem to contradict each other. 

This puzzle may now be understood as follows. The uncertainty relation ( [72D does 
not really follow from (|69D, because the time t, unlike x, is not an operator in quantum 
mechanics. However, according to (|7^) and the definition of ^-product, it is an effective 
description of the fact that the extent of spatial nonlocality (the shift in x for $ in 
the interaction term) is proportional to the energy scale, which is approximated by the 
free energy (—(9^ + m^)^/^ for small A. Incidentally, this relation is reminicent of the 



characteristic property of strings |2^. This is also one of the reasons why spacetime 
noncommutative field theories are of interest to us. On the other hand, in an exact 
treatment, the consistent high energy completion of spacetime noncommutative theory 
is presumably the NCOS theory, which contains a lot more degrees of freedom. 
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A One-Loop Unitarity for Derivative Interactions 

Consider a scalar field for example. In D dimensional spacetime, its propagator is 

A(x) = ^(xo)A+(x) + ^(-Xo)A-(x), (73) 
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where 



A±(x) = -^^^i^ J d''k5{e + m')e{±ko)e 



ik-x 



(74) 



Im 




Figure 2 

The left hand side of the diagram in Fig. 2 is 

t{xf-x,)=ig^F{df,df)F{dP,df^)A,{xf-x^^)/\,{xf-xf^)^^^^^^ (75) 

where the interaction is modified by a (Hermitian) function F of derivatives of two 
particles. For the special case of interactions in noncommutative field theories, 

F^d^^)^d^^)) = e¥^''^'^'^'"\ (76) 

By Fourier transform, 

S(p) = ig' j d^a:e-P-F(9)F(-9)Ai(x(i))A2(a;(2))L(.)_,.)=,. (77) 
If there is no time derivative in F, we find 

Im(S(p)) = yJ^^2 ip-h-k2)- 

■{6{po - uji{ki) - uj2{k2)) + 6{po + uji{ki) + uj2{k2))), (78) 

which is exactly the statement of unitarity. The same derivation can be extended to 
other one-loop diagrams, the point being that spatial coordinates do not have to be 
explicitly integrated over in the calculation. Thus we see that interaction terms with 
spatial derivatives preserve unitarity. The derivation breaks down when F contains time 
derivatives. The unitarity was found to be preserved for like-like noncommutativity in 
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